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ABSTRACT 

We present an analysis of probability distribution functions (pdfs) of column density in 
different zones of the star-forming region Perseus and its diffuse environment based on 
the map of dust opacity at 353 GHz available from the Planck archive. The pdf shape 
can be fitted by a combination of a lognormal function and an extended power-law tail 
at high densities, in zones centred at the molecular cloud Perseus. A linear combination 
of several lognormals fits very well the pdf in rings surrounding the cloud or in zones 
of its diffuse neighbourhood. The slope of the mean density scaling law (p)p oc L a 
is steep (a = —1.93) in the former case and rather shallow (a = —0.77 ± 0.11) in 
the rings delineated around the cloud. We interpret these findings as signatures of 
two distinct physical regimes: i) a gravoturbulent one which is characterized by nearly 
linear scaling of mass and practical lack of velocity scaling; and ii) a predominantly 
turbulent one which is best described by steep velocity scaling and by invariant for 
compressible turbulence (p)lu^/L, describing a scale-independent flux of the kinetic 
energy per unit volume through turbulent cascade. The gravoturbulent spatial domain 
can be identified with the molecular cloud Perseus while a relatively sharp transition 
to predominantly turbulent regime occurs in its vicinity. 

Key words: ISM: clouds - ISM: structure - submillimetre: ISM - turbulence - meth¬ 
ods: data analysis - methods: statistical 


1 INTRODUCTION 

On e of the tasks of the High Frequency Instrument (HFI; 
see lLamarre et al. IB) operating on the Planck satel¬ 
lite is to map the spatial distribution of Galactic molec¬ 
ular and diffuse gas. The data on dust emission at fre¬ 
quency of 353 GHz allow for tracing the general struc¬ 
ture of molecular clouds (MCs) and star-forming regions 
in a wide range of spatial scales - from vast zones that 
include MCs and their neighbourhood at scales of tens 
of parsec through their diffuse envelopes down to their 
dense and clumpy central zones. Such a study would con¬ 
tribute to the numerous attempts of different authors to 
model and explain the MC stru cture throu g h analysis of 
the column density d istribution (|Lpmbardi, Alves fe Ladal 

12010 ; iBallesteros-Paredes et alJl2012h . of discrete, hierarchi- 
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cal clumps 
approaches 


I Kauffmann et all l2010al fbl) or combining both 
(Beaumont et al.l 20121). 


The Perseus region is particularly interesting for the 
study of general MC structure due to several reasons: i) it 
is compact, apparently axisymmetric; ii) it does not over¬ 
lap with other fore- or background clouds or cloudlets, due 
to its high galactic latitude (6 ~ —20°); iii) it is associ¬ 
ated with several sites of recent or active star-formation; 
iv) the whole reg ion has been extensively investigated (see 
iBallv et al.ll20oT for a review). 


The shape of the column density probability distribu¬ 
tion function (hereafter, A-pdf) contains important infor¬ 
mation a bout the physics of star-forming regions. Dust ex¬ 
tinct ion dKainulainen et ail 120091 ; I LombardL _Alyes_jMAda 


2011) as well dust continuum studies (e.g. I Schneider et al. 


20131 ) demonstrate that the A-pdf in Galactic MCs is nearly 


lognormal in some cases or can be fitted through a combi¬ 
nation of lognormal and power-law functions, in other cases. 
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Figure 1. The selected probe zones, drawn on the Planck dust opacity map: (left) the ‘Perseus region’ (PR), containing Perseus MC 
(the central ellipse) and a high-extinction ‘core’; (right) zones of diffuse gas (labelled; see also Tabic [TJ) in the neighbourhood of the PR 
(unlabellcd) 


A p urely lognormal A-pdf is typical for turbulent medium 
(e.g. IVazauez-Semadenilfl994l ) as the standard deviation de- 
pends on the Mach number and the forcing of turbulence 
jFederrath. Klessen fe Schmidtl [20081 ) . On the other hand, 
numerical simulations of self-gravitating medium indicate 
that the development of a power-law (PL) tail is to be ex¬ 
pected at high d ensities and at timescales, comparable to 
the fr ee-fal l time (IKlesse nl200d: iKritsuhj^Jomian^^VVagnm 


l201ll : iFederrath fe Klessenl 120131 : ICirichidis et al.1 1201 4) . 
Therefore, the A-pdf can be used as an effective research 
tool to disentangle complex issues like scaling laws in turbu¬ 
lent MCs, energy equipartitions in star-forming regions and 
their relation to cloud evolutionary status. 

In this Paper we study and compare the behaviour of 
the A-pdf in embedded zones of Perseus region, containing 
Perseus MC, and in zones of its diffuse vicinity. The per¬ 
formed analysis allows for the reconstruction of the density 
scaling law and for the description of the physical regime in 
the Perseus region. 


Table 1. Statistics of the selected elliptical zones. 


Zone 

Nr. 

of pixels 

R 

[ P c ] 

{n)r 

{ 10 21 cm- 2 ] 

Mass 

[ 10 4 M 0 ] 

Embedded zones 

in the ‘Perseus region’ 


‘Core’ 

4663 

5.3 

10.1 

3.4 

1 (the MC) 

14498 

9.1 

7.4 

7.8 

2 

32299 

13.6 

5.2 

12.0 

3 

57339 

18.1 

3.9 

16.0 

4 

89649 

22.6 

3.1 

20.0 

5 

123812 

26.5 

2.8 

25.0 

6 (the PR) 

163351 

30.4 

2.5 

30.0 

Diffuse zones in 

the neighborhood 



dl 

155132 

30.4 

1.0 

14.7 

d2 

155121 

30.4 

0.8 

10.0 

d3 

147628 

30.4 

0.9 

10.0 

dl U d2 U d3 

318479 

41.4 

0.9 

19.6 


2 OBSERVATIONAL DATA AND PROBE 
ZONES 

We used the component map of dust opacity at frequency 
353 GHz (angular resolution: 1(72 per pixel) available from 
the Planck archivff] to extract a large area containing the 
Perseus MC. Dust emission data are more appropriate for 


1 http://irsa.ipac.caltech.edu/data/Planck/release_l/all-sky- 
maps/previews/HFI_CompMap_DustOpacity_2048_R1.10 


our study than dust extinction maps for two major rea¬ 
sons. First and foremost, extinction mapping depends on 
detection of background stars and thus suffers from selec¬ 
tion effects. Precise photometry of weak stellar sources is 
possible in zones of low to moderate extinction which biases 
extinction maps towards lower values. Consequently, masses 
of fragments of f ixed size can be essentia lly underestimated 
(cf. Sec. 2.2.2 in lKauffmann et al.ll2010bh . Second, Planck 
data offer the opportunity to probe regions of very high col¬ 
umn densities, not resolved in dust extinction maps due to 
a loss of background stars. 

Recent research indicates that a distance gradient exists 
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Figure 2. Column density pdfs in the PR (Fig. 1): (left) N is in absolute units, the lognormal fits of all pdfs (dashed) and the fit of the 
PL tail in the largest zone (solid) with its slope n are given; (right) N is normalized to its mean value in the considered zone; the span 
of is plotted with vertical shaded area. Only some of the pdfs are shown for clarity. 


across the Perseus region (ISchlaflv et aljfibldl i. Therefore all 
spatial sizes used in this Paper were calculated adopting a 
mean distance D = 260 ± 40 pc to the Perseus MC tak- 
ing into account th e parallax estimates to its western part 
llHirota et alJl201uh and th e variety of other ones to different 
parts of it (for review, see lBallv et al.l[2008l ). 

The selected zones to derive the 7V-pdf are shown in 
Fig. [T] while general physical information on them is pro¬ 
vided in Table |T] The zone labelled ‘Perseus region’ (here¬ 
after, PR) contains six embedded ellipses of increasing size 
and identical orientation angle, following the morphology of 
Perseus MC (Fig. [I] left). The central ellipse (#1) encom¬ 
passes the very cloud. The next two (#2 and #3) cover the 
cloud and also filamentary structures connected to it. The 
three largest elliptical zones (#4, #5 and # 6 ) include the 
cloud and its filamentary vicinity plus increasing areas of 
diffuse gas and of the south-eastern part of the Perseus OB2 
association, but without overlapping with the California MC 
(outside the northern edge of Zone # 6 ). The dense molecu¬ 
lar gas with N > 10 22 cm -2 is concentrated in Zone #1 (the 
Perseus MC) wherein we delineated also a ‘core’ with high 
dust extinction and as sociated with the Hu region G159.6- 
18.5 jBallv et al.ll2008l ). To compare the contribution of dif¬ 
fuse gas (atomic and molecular) to the IV-pdf, we selected 
additionally three zones in the western neighbourhood of PR 
with the same areas like the largest elliptical zone (# 6 , i.e. 
the whole PR) and spanning similar column density ranges 
(Fig. [U right; Tablc[T] bottom). 

A standard way to estimate column density from ob¬ 
servations of interstellar medium is from extinction mea¬ 
surements, adopting a constant ratio of column density 
N to visual extinction Ay. The value of N(Hi)/Ay 
which is fr equently used in different works is 1.9 x 
10 21 cm -2 (iBohlin et aljfl97gl ). It has been obtained from 
photometry of OB stars, with an upper limit in the chosen 


sample E(B — V) ~ 0T8 or, Ay < 2T5. In view of this limi¬ 
tation, transformations to column density in MC fragments 
of high extinction (typically populating the PL tail of the 
IV-pdf) should be considered with caution. Therefore, in this 
Paper, we derive N( H) directly from dust opacity T 353 at 353 
GHz, provided from Planck da ta. As evident from Fig. 6 in 
IPlanck Collaboration 19l (1201 il l, there is a very good linear 
correlation between dust opacity and total hydrogen column 
density N( H) obtained from Hi 21 cm maps, except for the 
range 0.7x 10 21 ^ V(H) ^ 5x 10 21 cm -2 wherein undetected 
(‘dark’) molecular gas may affect the results. Neglecting the 
latter effect, we adopted a linear conversion formula from 
dust opacity to hydrogen column density: 

N(¥L) = C 1 T 353 + Co , (1) 

where coefficients Co and Ci were obtained from fitting the 
IV(H)-binned representation of the correlation at column 
densities below 0.7 x 10 21 cm -2 (see Fig. 6 in Planck Collab¬ 
oration 19, 2011). The possible uncertainty of the calculated 
IV(H) is about a factor of 2 which is not significant for the 
analysis presented in this Paper. 


3 COLUMN DENSITY PDFS 
3.1 Perseus region 


The derived N-pdfs in the PR and in its embedded zones 
(Fig. [I] left) are plotted in Fig. [2] (left). Apparently, their 
shape can be described by a combination of a lognor¬ 
mal function around the peak and a power-law (PL) tail 
with slope varying in a narrow range between —1.7 and 
— 2.1. Such behaviour is in agree ment with other obser¬ 
vations of nearby MC complexes (Kainulainen et al.l 120091 : 
iFederrath fe KIessenll20l3 : ISchneider et al.ll2014al bh or with 
results from numerical simulations of self-gravitating clouds 
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dKritsuk. Norman fc Wagneijl20lH '). The slope and the lower 
column-density limit Ai, )w of the PL tail were calculated 
through the method Plfit dClauset et al.l 12009 1 which ex¬ 
tracts these characteristics from analysis of the unbinned 
observational data. In this way we avoid subjective estima¬ 
tions - note, that the obtained values N[' m (Tabic 0 are 
lower than one would put ‘by eye’ (say, at ~ 8 x 10 21 cm~ 2 ; 
Fig. [2] left). 

Normalisation to the mean column density (A) in each 
considered zone yields A-pdfs with practically identical 
shape and Nf >vr /{N) (Fig. [2] right), for a wide range of effec¬ 
tive radii. Physical interpretation of this result is suggested 
in Sect. U 


3.2 Diffuse zones 

The A-pdf in each diffuse zone exhibits a complex shape 
with several discernible peaks in whose vicinity the distri¬ 
bution is apparently lognormal (Fig. [3] left). The same be¬ 
haviour is found when one considers the A-pdf in any of 
the rings outlined by the boundaries of embedded elliptical 
zones in the PR, excluding the PL tail (Fig. [3] right). 


3.2.1 Method for pdf decomposition 

The column-density distribution in any of these regions of 
diffuse gas can be represented as sum of lognormal functions 
of type 






: exp | 


[igW ^)] 2 

2a? 


where the total number m of components (1 m) spans 

typically from a few to dozen and the parameters Ai, a; and 
a,i are obtained from the following fitting procedure of the 
whole distribution: 


• Prominent local peaks, i.e. of considerable width and 
height, are located and guess values A® are assigned to 
them. It is appropriate to start from such peaks with broad 
wings (large ai) since they have greater statistical weight. 

• The distribution is fitted in close vicinities of these 
peaks through single lognormals with fitting parameters 
N^\ and a- 1 - 1 . 

• The total fitting function YllLi lgn^ is composed. Then 
the whole distribution is fitted and next approximations of 
the fitting parameters N- 2 \ and a^ are obtained. 

• Small local peaks and inflexion points are located and 
the corresponding number of lognormal components lgn i is 
added to the fitting function, increasing m. The previous 
three steps are repeated and a next set n!> 3 \ a® and a!f' > 
is calculated for all components. 

• In case the achieved approximation of the distribution 
is not satisfactory, a few further components could be added. 
The latter have small statistical weights and contribute only 
for small local improvements of the fit. 

• Now the fitting-parameters’ space ( Ni,ai,ai ) (i = 
1 ,...,m) is severely restricted and the total fitting func¬ 
tion lg n i is fitted via NLLS Levenberg-Marquardt al¬ 
gorithm until the \ 2 criterion for fit’s goodness (at 95% con¬ 
fidence level) is satisfied. 


All obtained fits in the diffuse zones and in the rings 
of the PR are plotted in the Appendix m f^s. eu and lA 2 l 
The achieved fits’ goodness is given in Table [All 


3.2.2 Lognormal components and spatial scales 

A possible interpretation of this decomposition of the TV- 
PDF in a region to lognormals is that domains of various 
typical column densities Nj and velocity dispe rsions (related 
to ai; see e.g. Sect. 3.6 in lFederrath et aI1l20ld l contribute to 
it. Each domain is constituted of non-overlapping fragments 
of the considered region. Through the third parameter ai of 
a given component, we may define the effective size Li of a 
domain and hence introduce the notion of spatial scale. Each 
scale has its own lognormal column-density pdf lgn i (A) and 
fractional area f \gn i (N)/ JA j'lgn^TV). Then the effective 
sizes of scales included in a zone of diffuse gas with total 
area S and effective radius R = sjS/n are derived by: 



( 2 ) 


We stress that spatial scale is a rather abstract quantity 
which includes, in general, a set of fragments with various 
sizes (down to single pixels) and shapes. The notion is illus¬ 
trated in Fig. [4] which represents one possible spatial real¬ 
ization of four scales “detected” through decomposition of 
the A-pdf in an exemplary diffuse region. The larger scales 
L\ and L 3 (large ai/'^2 i m) consist of fragments spanning 
broad column density ranges (wide IV-pdfs) and have dif¬ 
ferent typical densities Ni. On the other hand, the smaller 
scales L 2 and L 4 (smaller ai/JA®*) are more homogeneous 
(narrow A-pdfs), with minor variations of A. 

The A-pdfs in the PR cannot be represented as com¬ 
binations of lognormals but note that their column density 
ranges almost entirely fall into the one, corresponding to the 
PL tail of the largest zone ff 6 (Fig. [3] left). Therefore one 
may consider zones ffl to as subscales of ff6, adopting 
their effective radii R as spatial scales. It will be shown be¬ 
low that all zones in the PR (except the cloud ‘core’) obey 
a mass scaling law with identical slope. 


4 PHYSICAL ANALYSIS 
4.1 Assumptions and definitions 

The analysis presented below is based on several reasonable 
physical assumptions. Some of them a re dis cussed and jus¬ 
tified in iDonkov, Veltchev fe Klessenl (l201lh , in the frame¬ 
work of a statistical approach for description of MC struc¬ 
ture at an early evolutionary stage. Here we briefly sum¬ 
marise them as follows: 

(i) Power scaling laws of velocity dispersion and me an den¬ 
sity (often called “first and second Larson’s laws”; iLarsonl 
ffl) , as well of mass: 

ul oc L^ , (3 > 0 . (3) 

(p)l oc L a , (4) 

MlocA. (5) 

The natural relation Ml = ( p)lL 3 leads to: 
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Figure 3. Examples of N -pdf (open symbols) decomposition to lognormals (dotted) in a diffuse zone (left) and in a ring (right) delineated 
by the boundaries of two successive elliptical zones in the PR (cf. Table [T}. In the latter case, the location of the lower boundary 
of the PL tail is indicated (dashed). The embedded graph demonstrates how the iV-pdf in a considered ring results from iV-pdfs of the 
elliptical zones (filled symbols). 


7 = a + 3. 


( 6 ) 


(ii) An equipartition relation between gravitational and tur¬ 
bulent energy of type E^ n ~ /|IT|, where / is a constant- 
The existence of such energy equipartition at different spa¬ 
tial scales and evolutionary stages in a MC is observationally 
verified and is explained by the coupling of various physical 
proc esses in the interstel l ar medium to each other (see Sect. 
1 in iHennebelle fe Falgaronel 2012j]_. It is supported as w ell 
by numerical simulations ( Vazquez-Semadeni et al1l2007l ). 

Hence one obtains a relation between the scaling indices 
of mass and velocity dispersion: 


1,, 2 3 

-M L u L oc -G— 


ILL OC 




i.e. jS = (7 — l )/2 


(7) 


How can one estimate the mean density at given spatial 
scale? The most straightforward approach is to define {p)l 
as a quotient of the mean column density and the scale. At 
a scale Li in the diffuse rings and zones, characterized by a 
lognormal A-pdf with peak Ni, (A ) is simply Ni. On the 
other hand, A-pdfs in zones of the PR are asymmetric and 
in this case we use the estimates of (A) from Table [T] (col. 
4) to calculate (p)jj « (N)/R. 

The mean-density scaling laws in the diffuse rings, in 
the elliptical zones and in the test simulation (see Sect. l4.3ll 
are derived from the A-pdfs via a weighted least-squares fit¬ 
ting described in Appendix [B] These scaling laws and those 
of velocity dispersion h old in the inertial range of a super - 
sonic turbulent cascade JPonkov, Veltchev fe KlesserfeOlU) . 
Therefore it is instructive to check the behaviour of the flux 
of kinetic energy per unit volume through the cascade in 
case of compressible turbulence, defined as: 


2 Note that the value of / is not important for this analysis, 
e.g. we do not assume gravitational boundedness (/ 5 0.5) or 
virialization of the considered volume. 



Figure 4. Schematic representation of spatial scales L, in our ap¬ 
proach. The border of a considered diffuse region is drawn (white 
line). The positions and shapes of the fragments are arbitrary and 
represent one possible spatial realization of the scales. 


/ _\ u i / „\ U L 

\P)l - = (p)L 77 — 

tl L/ul 



( 8 ) 


where tl is the turbulent crossing time at scale L. The pos¬ 
sible scaling of this quantity is part of the physical analysis 
presented below. 


4.2 Physical regimes in the Perseus region and its 
vicinity 

We suggest that the A-pdfs as well the underlying p-pdfs 
reflect two types of physical regime which govern the gen¬ 
eral structure of MCs. The PL tail corresponds to a regime 
in which gravity determines structure whereas the pdf part 
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Table 2. Parameters of the JV-pdfs derived in the embedded 
elliptical zones of the PR. Column densities are given in units 
[10 21 cm -2 ]. Notation: N m - peak of the distribution, n - slope 
of the PL tail, /pl(S) - PL area fraction, - PL mass 

fraction. 


Zone 

N max 

<JN 

A, pl 

low 

n 

/pl(S) 

/pl(M) 

‘Core’ 

5.97 

0.14 

5.95 

-1.67 

0.88 

0.84 

1 

4.51 

0.18 

5.22 

-1.84 

0.85 

0.77 

2 

3.35 

0.17 

3.90 

-1.83 

0.83 

0.72 

3 

2.15 

0.12 

2.85 

-1.85 

0.82 

0.71 

4 

1.81 

0.12 

2.39 

-1.82 

0.83 

0.71 

5 

1.80 

0.13 

1.98 

-1.93 

0.84 

0.73 

6 

1.89 

0.16 

1.85 

-2.09 

0.85 

0.74 


that can be decomposed to lognormal components corre¬ 
sponds to a regime with a prevailing role of turbulence. This 
concept is consistent with the fact that a well-developed 
PL tail, spanning about one order of magnitude or more, 
is present in pdfs of regions containing the very MC (Fig. 
[2j while a pdf with lognormal components and without PL 
tail is characteristic of the rings in the PR and of the diffuse 
zones in its vicinity (Figs. I All and IA2I) . The derived mean 
density scaling laws in the suggested two physical regimes 
also differ significantly - they are illustrated in Fig. [5] Let 
us examine the emerging picture closely. 


4-2.1 Zones forming the power-law tail of the N-pdf 

As already commented in Sect. 13.21 the embedded zones in 
the PR can be considered as subscales of the largest Zone ff 6 
and thus the latter should determine the balance of energies 
in the cloud structure. The A r -pdfs in these zones, normal¬ 
ized to the corresponding mean column densities (N)r, turn 
out to be strikingly identical (Fig. O right, and Table [2]), 
with the same width ai v, peak N i 11 s , x /(N)r and lower bound¬ 
ary AC/ (N)r of the PL tail. (Only the MC ‘Core’ exhibits 
an exception from this trend.) This result might be inter¬ 
preted by assuming that the Perseus MC is self-gravitating 
and its general structure is described by power-law surface- 
and volume-density profiles. 

Let a region with area S and volume V, centred 
at L = 0, is characterised by a surface-density profile 
S (L) = Eo(L/l pc)~ f and a volume-density profile p(L) = 
po(L/l pc)~ p . This yields scaling laws of column density 
(N ) oc L aN (ajv = —t) and mean density {p)l oc L a 
(a = — p), where those quantities are averaged over the en¬ 
tire profile. Then, if the column- and volume density pdfs in 
the region display PL tails, one gets for their slopes n and 
q, respectively: 


pjv(ln(E/Eo)) 
p(ln(p/p 0 )) = 


dS 

= d ln(E/E 0 ) 
dV 

d\n(p/p 0 ) * P 


oc E 


-3/p 


-2/t 


q = 


n = —2/t 
-3/p 


(9) 

( 10 ) 


_As demonstrated by iKritsuk, Norman fe Wagneil 

(1201 ll) . the relation between the volume-density and the 


surface-density profiles in isotropic medium is: 

+ oo 

E(L) = J p((L 2 + x 2 ) 1/2 ) dx = ... 

— OO 

+ oo 

oc L 1-p J (1 + (x/L) 2 )~ p/2 d(x/L) , ( 11 ) 

o 

where the integral is converging for 1 ^ p ^ oo. Obviously, 
t = p - 1 = —(a + 1). 

Now the derived mean density scaling index a = —1.93 
in the PR (Fig.0 top; blue symbols) leads through equation 
[9] to a PL-tail slope n = —2.15, in a very good agreement 
with the value —2.08 obtained directly from the observa¬ 
tional ALpdf (Fig. [5] left). On the other hand, one gets from 
equation [ 6 ] a mass-size scaling index 7 = 1.07. Recalling 
assumption (ii) for equipartition of energies, we derive a ve¬ 
locity scaling law of vanishing slope (cf. equation Q: 

u L oc L (l ~ 1)/2 oc L p , /3 « 0.04 (12) 

This surprising result is consistent with the dendro¬ 
gram analysis of cloud fragme nts from numerica l simula- 
tion s of self-gravitating c loud llShettv et al] l201(il . see Fig. 
2 in lstanchev et aLIl2013l i. Its physical meaning is that the 
sum of the specific kinetic and gravitational energies at each 
scale within the density span of the pdf tail is an invariant: 

1 2 3 Ml . , , 

-u L - - G— = mv(L) . (13) 

Note that this relation is not equivalent to the assumed 
equipartition of energies at each scale, although derived from 
it. The kinetic and gravitational energy terms contain as well 
the mass which is scale-dependent. 

From equation [13] one obtains for the flux of kinetic en¬ 
ergy (p)lu\/L oc L -3 . This finding means that the kinetic 
flux increases strongly at small scales which can be naturally 
explained with matter acceleratio n due to gravity. Ou r result 
is very different from the one of lGaltier fe Banerieel (120111 1 
who derived from theoretical considerations of compress¬ 
ible isothermal turbulence ( p)lu\/L = inv(L) for purely 
solenoidal forcing and ( p)lu\/L oc L 2 ^ 3 for strong compres- 
sive forcing. Th e latter relation was confirmed as well by 
iFederrathl (|2013l l from a number of high-resolution numeri¬ 
cal simulations of compressible turbulence without gravity. 
It is evident that the structure of the PR centred on Perseus 
MC cannot be modelled taking into account only interstellar 
turbulence; self-gravity is an essential factor in constructing 
the correct physical picture. 


4-2.2 Diffuse zones within the Perseus region 


The derived mean density scaling law in the rings of the PR 
is shown in Fig. [5] top (red symbols). It is shallower than the 
classical “second Larson’s law" with a = —1.1 for MCs and 
cloud fragments (lLarso Taking the density scaling 

index a ~ —0.80, one gets 7 = 2.20 and /3 « 0.60. Then - in 
comparison with the analysis in the previous Section, - we 
find here a scale-invariant flux of the total energy per unit 
volume: 


\\(p)lu l - 

L/ul 


inv(L) . 


(14) 
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Figure 5. Power-law scaling of mean density: (Bottom) in the 
diffuse zones outside the PR (small circles) and in their union 
(large circles); (Top) in the PR (filled circles) and in its diffuse 
rings (open circles). The slope in the former case was derived 
excluding the cloud ‘core’ (square). 


In particular, it follows that the flux of kinetic energy 
per unit volume is also invariant: 

3 

(p) L = inv(L) , (15) 

in agreement with simulations of interstellar isothermal 
turbulence w i thout gravity and with mixed forcing by 
iKritsuk et al.l (120071 1. 


4-2.3 Diffuse zones outside the Perseus region 

Mean density in the outer diffuse zones is found to scale with 
an index about the Larson’s value (Fig. [5] bottom). The con¬ 
sideration of their union dl U d2 U d3 yields the same result 
with a larger dispersion. We attribute the different scaling 
law in comparison with the PR rings to the actual variety of 
distances. All spatial scales and, respectively, mean densities 
in the diffuse zones outside the PR were calculated adopting 
D = 260 pc (to the Perseus MC) which could be far from 
the real values for a number of lognormal components, in 
view of the large sizes of the zones and the remoteness of 
some parts of them from the cloud. 


Table 3. Summary of the used simulation. 


Box size 

256 pc 

Simulation code 

FLASH 

Boundary conditions 

periodic 

Initial temperature 

T = 5 000 K 

Initial density 

1 cm -3 

Turbulence 

decaying 

Mach number 

0.4 

Magnetic field (aligned along the X axis) 

3 pG 

Refinement levels 

11 

Maximum resolution 

0.03 pc 

Jeans length 

8 grid cells 


4.3 Comparison with numerical simulations 

For a numerical test of the presented results, we ran a mag- 
netohydrodynamical simulation of MC evolution includin g 
gravity performed with a FLASH code dFrvxell et a! .Il20()7)h . 
The simulation properties are summarised in Table [3j The 
cloud formation was modelled through convergence of two 
large-scale cylindrical streams, each 112 pc long, with a ra¬ 
dius of 64 pc. They are given an initial supersonic inflow 
velocity (isothermal Mach number of 2) in a warm neutral, 
initially homogeneous medium and collide at the centre of 
the numerical box. Turbulence is not driven continuously 
and, at late evolutionary stages, is due to fluid motions, 
driven by gravity. Feedback by sink particles is not included 
in the simulation. 

Similar numerical set-up w as c hosen in the works o f 
IVazauez-Semadeni et al . ([200 71 and lBaneriee et ahl (l2009h . 
The former study showed the establishment of equipartition 
between gravitational and kinetic energy (cf. our physical as¬ 
sumption (ii) in Sect. IXT1) in dense gas domains (> 50 cm 3 ) 
and at evolutionary stages t > 10 Myr. The latter work re¬ 
vealed that such condensations form mainly in knots of in¬ 
tersecting filamentary structures and those located in the 
central region of the flow collision plane host sink particles. 
Now, if one rotates the flow collision plane from a line-of- 
sight angle 0 = 0° (face-on view) to 9 > 80°, it should 
be expected that a cloud similar to Perseus MC will be “ob¬ 
served”: an elongated, quasi-axisymmetric dense region with 
some filamentary ‘t runks’ in its near vic inity (cf. our Fig. [T] 
with Fig. 1, right in lBaneriee et al1l2009l h Indeed, this is the 
picture on the column-density maps obtained from our sim¬ 
ulation at evolutionary stages t > 15 Myr, plotted in Fig. [fj] 
left. 

4-3.1 Column-density pdfs 

We delineated a rectangular region of the same effective size 
like the PR (Zone #6 in Table [l]) and derived the IV-pdfs 
in it, varying the line-of-sight angle in the range 70 ^ 9 ^ 
88 °. Some of the obtained distributions are shown in Fig. 
[62 middle. The PL-tail range as detected by Plfit increases 
to about one order of magnitude at t > 18 Myr. The slope 
|n| gets shallower with evolution but stabilizes gradually at 
values 2.1 — 2.3 for 80 < 8 < 84°, in excellent agreement 
with our result for the PR (Fig. [5]). Below some critical value 
9 ~ 79° the PL tail becomes very steep and too short, i.e. 
AC increases discontinuously by a factor 2-3. On the other 
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Figure 6. Column-density pdfs (middle) and scaling laws of mean density (right) at three evolutionary stages of MC evolution (left). 
All column-density maps are produced for a fixed line-of-sight angle 9 = 84° and the considered rectangular region is indicated with 
solid white line. Shaded areas in the middle panel denote the location of the lower boundary of the PL-tail. The pdfs for critical 

values of 6 are plotted with dashed and dotted lines (see text). 


hand, nearly side-on views (9 > 85°) are characterised by 
much longer and shallower PL tails: |n| is close to unity. 

Considering only the non-PL part of the pdf (N < 
iV£ w ), we decomposed it to lognormal components through 
the procedure described in Sect. 13.2.11 All obtained fits in 
the selected region, at three different evolutionary stages, are 
plotted in the Appendix in i ig.tm while the achieved fits’ 
goodness is given in Table EH We derived the mean-density 
scaling law, adopting the definition of (p)l from Sect. 14711 
The results for the corresponding late evolutionary stages 
are illustrated in Fig. [6] right. 

Evidently, the mean density scaling index a ~ — 1 does 
not evolve within the considered time-frame as the uncer¬ 
tainty decreases. Following the approach from Sect. 14.2.21 
and 14.2.31 one obtains for the mass and velocity scaling in¬ 
dices 7 ~ 2.0 and /3 ~ 0.5, respectively. Hence, the total 
energy per unit volume turns out to be scale-invariant: 

\{p)lu l - = inv(L) , (16) 

while the kinetic energy flux per unit volume scales weekly 
{p)lu\/L oc L~ 1//2 . This behaviour is closer to the ex¬ 
pected one for com pressible turbulent medium with purely 
solenoidal forcing (iGaltier fo Baneried 1201 ll ) but demon¬ 


strates that gravity nevertheless affects gas dynamics. The 
result might be affected as well by the presence of magnetic 
fields in our simulation. 


4-3.2 Volume density pdfs 

The use of a numerical simulation allows for a direct study 
of the volume-density pdfs which correspond to those of 
column-density. An example of such p-pdf of the gas phase 
with densities > 5 cm -3 is shown in Fig. Q top. The Plfit 
procedure derives a very steep PL tail ( q ~ —3.7) spanning 
about an order of magnitude and starting close to the distri¬ 
bution maximum. However, the statistics at high densities is 
very poor: a few dozens out of 175000 pixels have densities 
over 800 cm' 3 , which are typical for dense regions in MCs. 
High-resol ution simulat i ons of self-g ravitating turbulent me¬ 
dia (e.g. ICollins et al.l l201ll : iKritsuk, Norman fe Wagner! 
l201ll) produce p-pdfs whose PL tail starts at least 1.5 or¬ 
ders of magnitude from the peak. Obviously, the PL tail of 
the p-pdf has not been resolved in our simulation. One is not 
able to derive the actual slope but the data hint at a value 
not far from q = —1.5 as shown in Fig. [3 top. Referring the 
reader again to the relations l9l and flOl note that q ~ —1.5 
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Figure 7. Pdf and scaling of volume density in the rectangular 
region delineated in Fig. [6] left. Top: Density distribution at t = 
20.2 Myr (6 = 84°) with the lower limit of the PL tail from the 
Plfit procedure (arrow) and the presumed actual one (dashed). 
The presumed slope of the PL tail is shown for reference. Bottom: 
Scaling law of mean density, averaged over three evolutionary 
stages. 


yields a PL-tail slope of the A-pdf n ~ — 2 and a = —2 
which is consistent with the results for the PR (Figs. [2] and 
[5] top) as well with the simulational A-pdfs (Fig. [ 6 ] middle). 

Ignoring the presumed PL-tail (at densities > 
600 cm -3 ), we applied the procedure of lognormal decom¬ 
position (Sect. 13.2.1 f) to the rest of the p-pdf. (The results 
are given in Appendix 0 ) Thus the notion of spatial scales, 
introduced in Sect. 13.2.21 was retained as well in the 3D 
case; with the necessary correction of the coefficient in equa¬ 
tion [2] to (a,i/ JY Oi) 1 ^ 3 . The peak value pi of each lognor¬ 
mal component was adopted as mean density of the cor¬ 
responding spatial scale. The number of lognormal compo¬ 
nents at each of the considered three evolutionary stages 
16.2 ^ t ^ 20.2 Myr is high enough to allow for search of 
a correlation between density and size. The derived mean 
density scaling law is shown in Fig. [7] bottom. The slope 
varies between ~ —1.0 and ~ —1.4 but with a large scat¬ 
ter. Nevertheless, the correlation is seemingly independent 
on the evolutionary time although the spatial scales change 
their effective sizes and mean densities. Therefore it is safe to 
assume and to derive a single, time-averaged density scal¬ 
ing law, taking into account all data. The obtained slope 
—1.33±0.26 is steeper but - in view of the large data scatter 
- in general agreement with the result from the simulational 
TV-pdf analysis. 


5 TWO DOMAINS IN THE PERSEUS 
REGION: A HYPOTHESIS 


Our analysis hints at the existence of different physical 
regimes in two spatial domains associated with the PR: 
gravoturbulent in the Perseus MC and its vicinity, described 
by equation 1131 and predominantly turbulent in the diffuse 
neighbourhood, described by equations [14] and 1151 Below 
we propose an idea to estimate the effective size L g t of the 
gravoturbulent domain and of the transition zone between 
the domains. 

Let L and C are the scaling variables in the gravotur¬ 
bulent and in the predominantly turbulent domain, respec¬ 
tively. For simplicity, the transition zone is considered as 
an extended shell encompassing the gravoturbulent domain, 
with effective size £ t r, he. with thickness dti = \J Lgt + £t r — 
Lgt. The velocity dispersion vc in this diffuse shell obeys the 
invariant flux of kinetic energy (equation 11511 and the mean 
density scales with index a ~ —0.80 (see Sect. 1^72721) . On 
the other hand, adopting 7 ~ 1 and j3 ~ 0 in the gravo¬ 
turbulent domain (see Sect. I4.2~TI) . one gets a linear scaling 
of mass Ml = Mo(L/l pc) and a lack of scaling of velocity 
ul = const(L) there. 

It is safe to assume that the time of consideration 
dt is very small in comparison to the turbulent crossing 
time T£ in the predominantly turbulent domain and, ac¬ 
cordingly, the mass i ncrease AMl of the gravotur bulent 
domain is negligible (IHennebelle fc Falgarong I2012I . Sect. 
2.4.1). A further reasonable approximation is that the mass 
inflow is stationary: dMc/dt = const(f) > 0 (Klessen & 
Hennebelle 2011). Then the derivative of the total energy 
Etot — ALlu^/2 — 3GMl /5L is: 

dEtot _ 1 2 AMl to Ml AMl 3~,Ml AL /175 

dt 2 L dt 5 L dt 5 L 2 dt 

where the first two terms in the right-hand side account 
straightforwardly for mass accretion and the last one - for 
the size change associated with the mass increase. 

In view of the negligible AMl, the main contributor to 
dEtot /dt is the derivative of the kinetic energy of the pre¬ 
dominantly turbulent domain: \dEtot/dt\ ~ ^v\\dMc/dt\. 
Now, replacing dL/dt = (X/M^AMl/ dt in equation 1171 and 
using |dMj,/dt| ~ \dMc/dt\, one obtains: 


1 2 


I 2 
\ 2 Ul 


|g^U| ( /-i)|g2^| 

5 LI I 5 1 pc I 


recalling our basic assumption (ii) for equipartition E^ n = 
f\W\ at each scale. From the derived mean density scaling 
law in the PR (Fig. [5] blue symbols) one calculates directly 
M 0 = M l = 1 P c = (p)l= 1 P c -pc 3 ~ 10 3 M 0 . Thus we arrive 
at an expression to estimate £ tr , assessing the coefficient 
/ in the gravoturbulent domain and the coefficient in the 
velocity scaling law vc = Vo(L/l pc)^ in the predominantly 
turbulent domain: 

T^~) ^ ; ft ~ 0 . 6 O (18) 

1 pc / 5 Vq 1 pc 

Typical values of / and Vo in star-forming regions 
are provided from the extensive recent study of the virial 
parameter a v j r in MCs and embedde d MC fragments by 
iKauffmann, Pillai fe Goldsmithl (l2013h . Those authors de¬ 
fine Qvir = g{p) 217km/11UI, where g(p) is a geometrical fac¬ 
tor depending on the density gradient with index p. They 
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Figure 8. Modelled sizes of the gravoturbulent domain L g t 
(solid line) and of the transition zone dn (grey area) as func¬ 
tions of /, fixing no = 0.8 km/s and /3 = 0.60 in the turbu¬ 
lent domain. The locations of the ‘Core’, Zone 1 and Zone 2 
in the limiting cases of / (minimal value - open symbols; maxi- 
mal value - fill ed symbols) are c alcula ted according to Fig. 1 in 
iKauffmann. Pillai fc Goldsmith! J2OI3I ): the bars show the exten¬ 
sion of the transition zones. See text. 


found 0.2 < a v i r < 2 and 0.7 < vo < 1.1 km/s for the 
majority of sampled objects with masses 10 3 to few times 
10 4 Mg which correspond (through Mo) to spatial scales 
1 — 30 pc considered in this Paper. For a variety of possible 
density profiles an appr opriate mean value of th e geometri¬ 
cal factor is g(p) ~ 1.4 dBertoldi fc McKeelll992 ) and hence 
0-07 < / < 0.7, depending on the mass of the MC fragment. 

Assuming that mass in the transition zone still scales 
linearly, like in the gravoturbulent domain, one gets an equa¬ 
tion for L g t: 


M(dtr) — ModtrOcA,[(Lgt + d t r) —-Lg t J, OL — —0.80 (19) 

The obtained modelled values of L g t and dtr as func¬ 
tions of / (or of the virial parameter), for a fixed vo = 
0.8 km/s, are shown in Fig. [5] Making use of the ob- 
servational correlation between virial parameter and mass 
(IKauffmann. Pillai fc Goldsmithl 120131 . Fig. 1) and of lin¬ 
ear scaling of mass Ml = Mq(L/1 pc), one can compare 
the prediction of the proposed model with the locations of 
the internal elliptical zones of the PR. Evidently, the gravo¬ 
turbulent domain is somewhat larger than Zone #1 which 
encompasses Perseus MC as a compact and apparently ax- 
isymmetric region. This domain determines the A-pdf char¬ 
acteristics (Fig. [2] right) and the density scaling law (Fig. [5]) 
in the PR considered as a whole, while these quantities in the 
set of rings are typical for diffuse, predominantly turbulent 
regions. It is interesting also to note that the column-density 
range A > 10 21 cm -2 in the gravoturbulent domain (see the 
pdf of Zone #T in Fig. [5] left) corresponds entirely to the 
molecular gas p hase in Perseus according to the study of 
I Lee et al.l (12012) . 


6 DISCUSSION 

6.1 Possible uncertainties of the density scaling 
law 

The derivation of density scaling laws could be sensitive 
to various factors. Here we examine briefly three of them: 
map resolution, uncertainties due to the distance gradient 
to the PR and uniqueness of A-pdf decomposition. Our re¬ 
sults were obtained by use of the original map resolution 
which corresponds to ~ 0.1 pc, i.e. about the size of dense 
clumps. We mimicked the effect of lower resolution using 
moving median smoothing of the column-density map, as 
described in Appendix [A] wherein the corresponding change 
in the pdf shape is illustrated (Fig. IA3I) . Although one or 
two additional components would be required in such a case 
to improve fit’s goodness, the derived slopes of the density 
scaling law (Fig. 0 practically do not change both in the 
embedded zones and in the diffuse rings. The distance un¬ 
certainties in th e whole PR can be a ssessed roughly since 
the estimates of lSchlaflv et al.l (120141 ) cover only the cloud 
itself and its close vicinity; Zones #T and #2. Linear fitting 
of those data yields an approximately constant distance gra¬ 
dient whose vector is nearly collinear to the large axis of the 
elliptical regions. Therefore one may consider the whole re¬ 
gion as being rotated at an approximately constant angle ip 
around the small axis. Such geometry would affect the ef¬ 
fective radii R and sizes Li (see eq. 0 by factor cos ip and, 
respectively, would shift the zero point of the derived den¬ 
sity scaling laws (Fig. 0 but would not affect the scaling 
indices. 

Next, one might question the uniqueness of A-pdf de¬ 
composition to a small number of lognormals. Is it possible 
to fit the observational column-density distribution by sum 
of another set of lognormals which shows up different phys¬ 
ical features? Such arbitrariness of the suggested approach 
is severely limited through its starting steps: 1) setting the 
number of components equal to the number of prominent 
local peaks of the A-pdf and, 2) choosing such guess val¬ 
ues of (Ai, a;, a,i) which nearly reproduce the location of 
each peak and the A-pdf shape in its vicinity. Finally a 
simultaneous Levenberg-Marquardt optimisation of all 3A/” 
parameters is performed to obtain the fitting curve. On the 
other hand, the termination of the fitting procedure depends 
on estimation of the fit goodness - in principle, the latter 
can be improved by adding one or two further components 
to account for single tiny fluctuations of the distribution. 
However, such additional component(s) would be of small 
effective size (small ai/T,idi) in comparison with the bulk 
of components and could slightly affect the scaling of mean 
density (cf. Fig. 0. The defined statistical weights (equa¬ 
tion IBII) also exclude essential changes of the mean density 
scaling index with adding up small components. 


6.2 Density scaling laws and physical regimes 

It may seem surprising that the derived mean-density scal¬ 
ing law in the diffuse rings a ~ —0.80 differs substantially 
from the one in the PR whereas - expressed in terms of 
mass-size relationship (7 ~ 2.2, equation 0, - it resem¬ 
bles that of a large sample of Galactic MCs with sizes from 
a few to several tens pc studied by iRoman-Duval et al.l 
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d2010l 'l. In fact, such behaviour is to be expected if one 
differentiates internal general structure of individual clouds 
from consideration of a sample of clouds. Total masses of 
MCs correlate well with th eir sizes following a single power 
law with 2 ~ 7 < 2.4 llLombardi. Alves fe Ladal l20ld : 
iRoman-Duval et al.M2010l ) while the index 7 within a cloud 
increases reversely to the scale L. The latt e r phe nomenon 
was discovered bv lLombardi, Alves fc Ladal (|2010| . see Fig. 
2 there) from a dust extinction study of 11 Galactic MCs and 
explained within a statistical mode l of general MC structure 
JPonkov, Veltchev fe Klessenll201lh . assuming equipartition 
of energies at each scale, characterized by its own pdf. The 
mean-density scaling law in the PR derived in this Paper 
confirms those results extending their applicability to the 
vicinity of the cloud itself (Fig. [5] L > 10 pc) while the 
value of 7 in the dense cloud regions (within the ‘Core’) 
might increase. A further investigation by use of Planck 
data, including other MCs with different general structure 
and various physical conditions (e.g. star-forming activity), 
would shed light on this important issue. 


On the other hand, the diffuse rings evidently obey 
a mass scaling law like that for single MCs. The power- 
law index 7 ~ 2.2 is also in agreement, within the 

uncertainties, with t he measurements from simulations 
jKritsuk_e£_aL 2007|; iFederrath, Klessen fe Schmidtl 120091 : 
iKritsuk et al .11 2001)1 ). This result, combined with our assump¬ 
tion (ii) for equipartition between gravitational and kinetic 
energies, lends support to the hypothesis that supersonic 
compressible turbulence dominates the diffuse gaseous struc¬ 
tures in cold interstellar medium as described by equations 
M and m It should be pointed out that assumption (ii) 
refers to abstract scales as defined in Sect. l3.2l and not to con¬ 
nected regions with well-defined boundaries for which sur¬ 
face energy terms are comparable to the volume ones. Never¬ 
theless, it is noteworthy that the mean-density scaling law of 
abstractly defined scales - both from observational data on 
the PR and from the numerical simulations (Sect. m - are 
close, within uncertainties, to a ~ —1.1 derived by Larson) 
(Il98lh for MCs and cloud fragments delineated as connected 
regions. Moreover, the striking coincidence between the de¬ 
rived scaling indices from column- and volume-density pdfs 
from our simulation confirms this result and supports the 
assumption of isotropy in turbulent gas medium. Thus the 
so defined spatial scales reproduce the fractal MC structure 
as derived by use of classical clump-finding techniques. 


The mean-density scaling law a ~ — 1 obtained from 
the performed simulations is slightly steeper than the one 
in the diffuse rings in the PR and seems to indicate slightly 
different physical conditions. One may attempt to integrate 
the obtained results for the PR and those from the simula¬ 
tions in a single scenario, interpreting the variety of physical 
invariants and scaling laws as referring to different evolution¬ 
ary stages. A clear indication for this is the behaviour of the 
flux of kinetic energy per unit volume ( p)lu\/L (equation 
The latter quantity is scale-independent in the diffuse 
rings of the PR, scales as L _1,/2 in the diffuse components 
of the simulational pdfs and scales as strongly as L~ 3 in the 
elliptical zones of the PR. That points to the increasing role 
of gravity in the gas dynamics. Possibly, the assumed three 
physical regimes hint at different self-similar gas structure. 


6.3 Relation to the velocity scaling law 


The density scaling laws in the PR are derived only by use 
of the procedure for A-pdf decomposition, without any ref¬ 
erence to the velocity scaling law. The relation between the 
two scaling laws relies essentially on our second physical 
assumptions for equipartition between gravitational and ki¬ 
netic energy (Sect. 14.11) . Introducing the latter, we aimed 
to explore the physics which governs the column-density 
statistics in star-forming regions and not to derive velocity- 
size relationships out of column density data. Nevertheless, 
a comparison of the (indirectly) estimated scaling index /3 
with observational and numerical studies might be instruc¬ 
tive. 


The l inewidth-size relat ions in the Perseus MC, ob¬ 
tained bv IShettv et alj (12012T ) from 13 CO COMPLETE sur¬ 
vey data, are plotted in Fig. [9] The considered objects are 
hierarchica l clumps, extracted by use of the Dendrogram 
algorithm (IRosolowskv et al]|2008l 'l. Therefore they are very 
appropriate for comparison to our embedded elliptical zones. 
Unfortunately, their size range overlaps with our study only 
within Zone #T and its ‘Core’ (cf. Fig. 0 up) so that we 
are able to present only a qualitative analysis here. A den¬ 
drogram of clumps is a multilevel hierarchy in which large 
embedded structures as well si n gle cond ensations are distin¬ 
guished (see Fig. 1 in lShettv et al.ll2012l 'l. It has been discov¬ 
ered from analysis of simulations of self-gravitating turbu¬ 
lent medium that series of massive dendrogram structures 
are characteri zed by unique mass-d ensity and velocity-size 
relationships (IStanchev et al.l 120131 1 . Such structures and 
their unique velocity-size relationships are shown in Fig. [5] 
with short dotted lines. They are representative for scales 
L > 2 pc within the Perseus MC. The averaged scaling index 
(/3) is in a rough agreement with the estimation from results 
in this Paper for the ‘Core’ and Zone #1. A tendency of 
shallowing in regard to the embedded zones outside Perseus 
MC is sensible, although the lack of molecular-line data for 
larger scales leaves the issue open. 

The simulation used in this Paper to compare with 
the derived density-scaling laws is not relevant to test the 
velocity-size relationship. It implements decaying initial tur¬ 
bulence which cannot provide a large inertial range to de¬ 
rive the velocity scaling law, although accretion on proto- 
stellar objects might be an effective driving mechanism at 
small scales. Nevertheless, the initial turbulence shapes sig¬ 
nificantly the density field in earlier evolutionary stages and 
hence is crucial for the cloud morphology. 

The virtual lack of velocity scaling in the PR (equation 
o indicates the essential role of gravity. It is noteworthy 
that the estimated effective size of the gravoturbulent do¬ 
main is about the size of the Perseus MC (cf. Table [T] and 
Fig- ID) f° r a variety of conditions: from weak to strong grav¬ 
itational boundedness (0.25 < / < 0.5). The size dti of the 
transition zone between gravoturbulent and predominantly 
turbulent regimes depends mainly on the balance between 
.Ekin an( i |1U| in the former regime and of coefficient vo of 
the velocity scaling law in the latter one. The typical values 
of d tr are about 10 % of size of the gravoturbulent domain, 
irrespective of the degree of gravitational boundedness. 
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Figure 9. Velocity-size diagram of large em bedded (dendrogram ) 
structures in the Perseus region studied bv IShettv et al.l 12012?) . 
with their unique scaling laws (dotted). The index 3 in the range 
of overlapping with the present Paper is estimated through the 
density scaling index a and shown for comparison. 
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6.4 The role of magnetic fields 

Finally, we briefly consider the possible effect of magnetic 
fields which were not taken into account in the adopted phys¬ 
ical frame (Sect. 14.111 . Detailed comparison between cloud 
evolution from simulations in non-magnetized and in moder¬ 
ately magnetized medium shows that the magnetic support 
delays the gravitational collapse whereas the corresponding 
stages of the iV-pdf evolution are identical, differing only by 
a small time factor dHennebelle et aPl2008l ) . About the time 
when gravitational collapse of the whole cloud is to start 
(i.e. when equipartition between gravity and turbulence is 
achieved), the shape of the column-density pdf in both cases 
is a combination of lognormal function pl us power-law tail 
at N > few times 10 21 cm -2 (see Fig. 1 in iHennebelle et al.l 
120081 ). Therefore one could not expect that the presence of 
magnetic fields would affect generally the presented physical 
analysis, given than the exact evolutionary stage of the PR 
is not known. 


7 CONCLUSIONS 

In this Paper, we presented an analysis of the column- 
density pdfs in the star-forming region Perseus and its dif¬ 
fuse environment using the dust opacity map at 353 GHz 
available from the Planck archive. Due to its high galac¬ 
tic latitude and approximately axisymmetric form, free of 
contaminating fore- or background cloudlets, the region is 
an appropriate test case for studying the general structure 
of molecular clouds and the physical processes that govern 
their evolution. On the other hand, the Planck data on 
dust emission offer the opportunity to probe regions of very 
high extinction, not resolved on dust extinction maps. 

Our main conclusions are as follows: 

1. The pdf shape can be fitted: i) by a combination of a 
lognormal function and an extended power-law tail at high 
densities, in zones centred at the molecular cloud Perseus; 
and, ii) by a linear combination of several lognormals, in 
rings surrounding the cloud or in zones of its diffuse neigh¬ 
bourhood. In the first case, the power-law tail of the pdf 
in the largest zone includes the pdf density ranges of all 
others (embedded) zones whose sizes may be considered as 


subscales of this zone. In the second case, each lognormal 
component of a given pdf is interpreted as a contribution 
of a separate spatial scale in the turbulent cascade. In that 
way, the notion of spatial scales L is introduced which al¬ 
lows for derivation of power-law scaling laws of mean density 
( p)l oc L a and mass Ml oc L 7 in the studied regions. 

2. The derived scaling laws in zones centred on the 
Perseus MC and those in zones of its diffuse vicinity or neigh¬ 
bourhood differ substantially: a = —1.93 (7 = 1.07) in the 
former case and a = —0.77 ± 0.11 (7 = 2.23) in the latter 
case. Assuming an equipartition relation between gravita¬ 
tional and kinetic (turbulent) energy at each spatial scale 
and a power-law scaling of velocity dispersion ul oc L 13 , this 
result bears evidence of two distinct physical regimes in the 
Perseus region: 

- gravoturbulent, characterized by nearly linear scaling 
of mass and practical lack of velocity scaling. 

- predominantly turbulent, characterized by a steep ve¬ 
locity scaling (/? ~ 0.60) and by invariant for compressible 
turbulence ( p)lu\/L , describing a scale-independent flux 
of the kinetic energy through turbulent cascade. 

General identification of these regimes with spatial domains 
in the Perseus region is physically sustained in point 5 bel¬ 
low. 

3. The mean-density scaling law a ~ — 1, derived both 
from column- and (within the a large data scatter) volume- 
density pdfs from the performed magnetohydrodynamical 
simulations with gravity, i s in e xcellent agreement with the 
classical study of iLarsonl dl98lh . This is striking since the 
result was obtained for abstractly defined spatial scales - 
not for connected regions in MCs, delineated by use of some 
clump-finding technique. 

4. The obtained three different physical invariants (resp., 
scalings of (p)lu\/L) in the PR, its diffuse rings and in the 
simulated cloud region seem to point to three different types 
of self-similar structure. 

5. Modelling the gravoturbulent spatial domain as cen¬ 
trally symmetric, with effective size L gt , and the predomi¬ 
nantly turbulent domain as its extended shell, we estimate 
Lgt and the thickness of the transition zone dt r - For rea¬ 
sonable values of the virial parameter, one obtains L g t ~ 
10 — 11 pc and dtr ~ 0.6 — 1.5 pc, i.e. the gravoturbulent 
domain can be generally identified with the Perseus MC in 
whose vicinity a relatively sharp transition occurs to pre¬ 
dominantly turbulent physical regime. The latter roughly 
coincid es with the tra nsition from molecular to atomic gas 
phase dLee et al ■ 20l3). 
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Figure Al. Decomposition of the TV-pdf (open symbols) in the 
diffuse zones to lognormals (dotted). The fitting function is plot¬ 
ted with solid line. 


APPENDIX A: DECOMPOSITION OF THE 
PDF: LOGNORMAL COMPONENTS AND FIT 
GOODNESS 

Below we show the pdf decomposition to lognormal compo¬ 
nents in the diffuse zones (Fig. I All) and in the rings out¬ 
lined by the boundaries of embedded elliptical zones in the 
Perseus region (Fig. IA2|) . For comparison, the pdf decompo¬ 
sitions in the simulated region at three different evolutionary 
stages (cf. Fig. [6]) are plotted in Fig. IA4I The goodness of the 
obtained fits of observational and numerical pdfs are given 
in Tables m and IA2I respectively. We studied the possible 
effect of insufficient resolution on the pdf, applying the tech¬ 
nique of moving median smoothing on the original Planck 
map: the value of N in each pixel was replaced by the median 
value in a square containing the neighbor pixels. A typical re¬ 
sult is shown in Fig. lA3l As expected, a lower map resolution 
would strengthen local peaks in N (i.e. the densities with 
higher probability) - the pdf becomes less smoothed and 
could require one or two additional components to achieve 
a satisfactory fit goodness. However, this additional compo¬ 
nents do not influence practically the derived density scaling 
law (slope changes of a few hundredths). 

The decompositions of the volume-density pdf in the 
simulated region at three different evolutionary stages are 
plotted in Fig. IA5I and the goodness of the corresponding 
fits is given in Table lA3l 

Evidently, the goodness of the obtained fits is satisfac¬ 
tory - in most cases (except for Ring 3-2 in the PR. and 
at t = 16.24 Myr in the simulations), the x 2 value is con¬ 
sistent with the two-sided test at 95% confidence level. As 
one can expect, the PL tail in simulational X-pdfs grows in 
terms of column-density range and statistical weight with 
the MC evolution (Fig. IA4l) . In contrast, the density-scaling 
law derived from the non-PL part does not change, with 
slope a ~ — 1 and uncertainty which decreases distinctly 
with evolution. 


Table Al. Goodness of the obtained fits to N-pdfs in the rings 
of the PR and in the diffuse zones. Notation: .Y|,j n - number of 
bins, m — total number of lognormal components, DF = degrees of 
freedom (Xb; n — 3m), x 2 rit low ( up ) “ lower (upper) critical value of 
the y 2 -distribution at 95% confidence level. The x 2 values, which 
are consistent with the two-sided test Xcrit.low ^ X 2 < X 2 r it,u P > 
are put in bold. 


Region 

Xbin 

m 

DF 

x 2 

Y 2 . 

/V crit, low 

-^crit, up 

Rings in the PR 






2-1 

19 

4 

7 

2.99 

2.17 

14.07 

3-2 

20 

2 

14 

1.58 

6.57 

23.69 

4-3 

22 

4 

10 

16.91 

3.94 

18.31 

5-4 

22 

4 

10 

4.68 

3.94 

18.31 

6-5 

26 

4 

14 

8.96 

6.57 

23.69 

Diffuse zones in the neighborhood 




dl 

55 

5 

40 

42.53 

26.51 

55.76 

d2 

51 

6 

33 

34.96 

20.87 

47.40 

d3 

52 

5 

37 

51.17 

24.08 

52.19 

dl U d2 U d3 60 

7 

39 

51.80 

25.70 

54.57 


Table A2. Goodness of the obtained fits to N-pdfs at three evolu¬ 
tionary stages of simulated MC evolution. The line-of-sight angle 
is denoted by 8. Other notation is the same like in Table lATl 


8 

Xbm 

m 

DF 

x 2 

•^•crit, low 

^•crit, up 




t = 

16.24 Myr 



82 

58 

12 

22 

31.70 

12.34 

33.92 

84 

54 

12 

18 

22.98 

9.39 

28.87 

85 

54 

13 

15 

21.47 

7.26 

25.00 

87 

55 

12 

19 

26.06 

10.12 

30.14 




t = 

18.59 Myr 



80 

31 

8 

7 

1.55 

2.17 

14.07 

81 

32 

11 

11 

12.53 

4.58 

19.68 

82 

32 

8 

8 

4.42 

2.73 

15.51 

84 

33 

8 

9 

10.71 

3.33 

16.92 




t = 

20.19 Myr 



81 

33 

10 

3 

1.32 

0.35 

7.82 

83 

34 

11 

1 

0.98 

0.01 

3.84 

84 

34 

11 

1 

1.57 

0.01 

3.84 

85 

35 

11 

2 

5.93 

0.10 

5.99 


APPENDIX B: DERIVATION OF DENSITY 
SCALING LAW THROUGH A WEIGHTED FIT 

What should be an appropriate and physically reliable mea¬ 
sure for the statistical weight Wi of each scale Li? Refer¬ 
ring the reader to the suggested procedure (Sect. 1372.11) , we 
point out that larger components have major contributions 
to the pdf fit while small components lead to minor local 
improvements. These contributions are proportional to the 
parameter a t and must be taken into account as one derives 
the scaling law of ( p). Hence, the chosen estimate of the 
statistical weight of a scale is: 
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Table A3. Goodness of the obtained fits to volume density pdfs 
at three evolutionary stages of simulated MC evolution and line- 
of-sight angle 0 = 84°. The notation is the same like in Table 

[ATI 


t[Myr] N hin m DF X 2 X«it, low *crit, up 

16.24 45 11 12 2.93 5.22 21.03 

18.59 42 10 12 6.57 5.23 21.03 

20.19 42 7 21 11.84 11.59 32.67 


Ay Ay 



N(x 353 )[10 21 cm' 2 ] 


Figure A2. Decomposition of the TV-pdf (open symbols) to log- 
normals in rings delineated by the boundaries of two successive 
elliptical zones in the PR (cf. Table [T}. The locations of the lower 
boundaries TVj^ G f the PL tail are indicated (dashed). Other no¬ 
tation is the same like in Fig. IaTI 
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t=20.19 Myr 
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Figure A4. Decomposition of the TV-pdf (open symbols) to log- 
normals in the simulated regions at three evolutionary stages of 
MC evolution (cf. Table [A~2l ). 


Figure A3. Illustration of the effect of map smoothing (mimick¬ 
ing lower resolution) on the pdf decomposition. See text. 


Wi(Li) oc a,i /^ at oc a/ y>. (bi) 

i i 

Examples of fits derived in the diffuse rings of the 
Perseus region and in the simulated cloud are shown in Fig. 
IBTI Evidently, the adopted weighting yields similar slopes 
to those calculated from standard NLLS procedure. 
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Figure A5. Decomposition of the volume-density pdf (open sym¬ 
bols) to lognormals in the simulated regions at three evolutionary 
stages of MC evolution (cf. Table ED . 
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Figure Bl. Correlations between effective size and mean density 
in linear (left) and log-log plot (right). Fits derived by standard 
NLLS procedure (dotted line) and by adopted weighting (solid 
line) are shown. 








